A0-A165  4*6 
UNCLASSIFIED 


ASVHPTOTIC  ANALYSIS  OF  A  TURBULENT  BOUNDARY  LAVER  IN  A 
STRONG  ADVERSE  PRESSURE  GRADIENT (U>  DAYTON  UNIV  OH 
RESEARCH  INST  H  B  BUSH  ET  AL.  JUL  87  AF0SR-TR-87-B9G2 
F49428-85-C-0137  F/G  2B/4 


1/1 


NL 


AD-A185  406 


1*.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


UHS  Bl£ 


\JSJ  i 


REPORT  DOCUMENTATION  PAGE 


2a.  SECURITY  CLASSIFICATION  g-  CT  E  | 

2b-  OtCLASSIFICATION/OOWNGlfc^O^|AyiQ.E|  j 


4.  PERFORMING  ORGANIZATION 

UDR-TR-87-35 


nr  NUMBER(S) 

^OD 


1b.  RESTRICTIVE  MARKINGS 

5' 3.  OfSTRI  BUT  I  ON/A  VAI  LABILITY  OF  REPORT 

_  Approved  for  Public  Release; 
Distribution  is  Unlimited. 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER'S) 

AFPSR-TR-  8  7-  0  962 

ISO L  7«.  NAME  OF  MONITORING  ORGANIZATION 

ftFOSR I  urn 


6.  NAME  OF  PERFORMING  ORGANIZATION  3b.  OFFICE  SYMBOL  .7*.  NAME  OF  MONITORING  ORGANIZATION 

University  *  U.yton  AFOCft  | 

•c.  AOORESS  (City.  Slat*  and  7.1?  Cod*)  7b.  ADDRESS  (City.  Stall  and  Zl?  Cod*) 

°f  Dayt°n  AFOSR/NM^bl 

300  College  Park  Ave  Bolling AFBDC  20332-6448 

Dayton,  OH  45469-0001 

to.  NAME  OF  FUNOING/SPONSORlNG  8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  lOENTIFICi 

ORGANIZATION  force  at  appUeabt*) 

Office  of  Scientific  Research _ NM _ F  49620-85-C-0137 _ 

to  AOORESS  (City.  Slat*  and  ZIP  Cod*)  10.  SOURCE  OF  Funding  NOS. 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

NM _ F  49620-85-00137 _ 

10.  SOURCE  OF  FUNDING  NOS. 


PROJECT 

NO. 


8c.  ^DORESS  (City,  State  and  ZIP  Code)  10.  SOURCE  OF  FUNDING  NQS.  _ 

J^vlX  \  PROGRAM  .  FRO-ECT  TASK 

Bolling  AFB,  DC  20332  elementno.  no.  no. 

11.  TITLE  (Include  Security  Clout  fleet  tout  Asymptotic  Analysis  IbiioaF  a504  fl-3 
of  a  Turbulent  Boundary  Layer  in  a  Strong _ _ 


WORK  UNIT 
NO. 


12.  PERSONAL  AUTHOR(S) 

W.B.  Bush,  and  L.  Krishnamurthv 

13.  TYPE  OF  REPORT  13b.  TIME  COVERED 

Interim  from  1  Jan  §7- 

19.  SUPPLEMENTARY  NOTATION 


14.  DATE  OF  REPORT  (Yt  .  Mo..  Do»  15.  PAGE  COUNT 

1987,  July  55 


fl7.  COSATI  COOES 

18.  SUBJECT  TERMS  ^Continue  on  rwverte  if  nece story  and  identify  by  block  number !  j 

FIELD  GROUP 

SUB.  GR. 

Adverse  Pressure  Gradient 

Mixing-Length  Closure 

Asymptotic  Analysis 

Reynolds-Averaged  Equations 

\ 

Limit-Process  Expansions 

Turbulent  Boundary  Laver 

19.  ABSTRACT  t  Continue  on  reverse  if  necettary  and  identify  by  block  number  / 

The  structure  of  an  incompressible  turbulent  boundary  layer  subjected  to  a 
strong  adverse  pressure  gradient  is  studied  by  means  of  an  asymptotic  analysis 
of  the  Reynolds  time-averaged  equations.  Limit-process  expansions  developed  in 
the  limit  of  large  Reynolds  number  reveal  a  relatively  thick  nondefect  layer  in 
the  outer  region  of  the  boundary  layer  near  the  exterior  inviscid  flow,  and  a 
relatively  thin  layer  near  the  wall.  To  leading  orders  of  approximation,  the 
momentum  balance  involves  convection,  pressure  gradient,  and  turbulent  stress  in 
the  outer  layer,  and  pressure  gradient,  and  turbulent  and  viscous  stresses  in  — ^ , ..  .7 

20.  DIST  Rl  BUTION/AV  AILABI  LIT  Y  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFI60/UNUMITE0  53  SAME  AS  RPT.  □  OTIC  USERS  □  UNCLASSIFIED 

22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  NUMBER  22c  OFFICE  SYMBOL 

< Include  Area  Codet 

Major  John  Thomas,  Jr. _ (202)  767-5026  AFOSR/NM 

DD  FORM  1473,  83  APR  EDITION  OF  1  JAN  73  IS  OBSOLETE. 


22c.  OFFICE  SYMBOL 

AFOSR/NM 

SECURITY  CLASSIFICATION  of  this  page 

r\  n (A 


97  £  H  o&o 


11.  (cont'd)  Adverse  Pressure  Gradient 

19'7'  the  inner  layer.  The  asymptotic  expansions  for  these  two  layers  are  matched  in 
an  arbitrary  intermediate  region,  wherein  the  streamwise  velocity  has  a  square- 
root  dependence  and  the  Reynolds  stress  has  a  corresponding  linear  dependence  on 
the  normal  coordinate.  The  leading-order  approximations  for  the  outer  and  inner 
layers  give  rise  to  similarity  formulations?  from  which  appropriate  similarity 
formulations  for  the  distinguished  intermediate  layer  have  been  identified  and 
developed.  These  latter  formulations  are  employed  to  analyze  available 
experimental  data.. 


SSCUAITV  CLASSIFICATION  OF  THIS  F AGl 


ASYMPTOTIC  ANALYSIS  OF  A  TURBULENT  BOUNDARY  LAYER 


E 

r 

E 


IN  A  STRONG  ADVERSE  PRESSURE  GRADIENT 


W.  B.  Bush  and  L.  Krlshnaaurthy 


E 


C 


JULY  1987 


4 


E 

L 

L 

r 

L  UNIVERSITY  OF  DAYTON 

RESEARCH  INSTITUTE 
DAYTON,  OHIO  45469-0001 

L 


L 

L 


AFOSR.TR.  37 


0  962 


PREFACE 


The  completion  of  this  research  and  the  documentation  of  the  results 
were  supported  by  the  Research  Funds  of  the  Office  of  Director,  University  of 
Dayton  Research  Institute.  Dr.  William  B.  Bush,  Senior  Associate,  King.  Buck  & 
Associates,  Inc.,  2384  San  Diego  Ave.,  Suite  2,  San  Diego,  CA.  92110,  and  Dr.  L. 
Krishnamurthy,  Senior  Research  Engineer,  Applied  Physics  Division,  University  of 
Dayton  Research  Institute  express  their  appreciation  to  Drs.  John  Murat  and 
Eugene  Gerber  for  their  encouragement.  One  of  us  (L.  K.)  was  supported  in  part 
by  the  U.  S.  Air  Force  Office  of  Scientific  Research  under  Contract  No.  F-49620- 
85-C-0137  (Major  John  Thomas,  Jr.,  Project  Monitor).  The  authors  are  Indebted 
to  Ms.  Ellen  Bordewisch  for  her  skillful  assistance  in  the  preparation  of  this 
report.  The  authors  are  pleaned  to  note  that  the  present  analysis  and  results 
have  benefited  from  earlier  interactions  with  Drs.  Noor  Afzal  and  Francis  E. 
Fendell. 


% 


'6 ; 


'j .  'c 


°a  % 


Accesio 

For  T  * 

NTIS 

cut  :<l  \J  ! 

D1!C 

ta-<  n  i 

Uiijr.:  c 

. «  •  '  1 

.  .  J 

J(»  -tit 

i 

By 

i 

L«: in 

1*1  ■  .  1 

•’  1  •  '  ^ 

,  ■■■-.-—•  1 

n.  ♦  * 

^  | 

A  /' 

!  1  I 

/)-/ 

.  i 

:  i 

r ' } 


IrWv 


1.  INTRODUCTION 


The  vast  literature  that  addresses  turbulent  boundary  layers  developing 
under  significant  longitudinal  (favorable  or  adverse)  pressure  gradients  attests 
to  their  (practical  and  f undanental )  iaportance.  Much  is  known  about  these 
flows  from  the  experlaental .  theoretical,  and  computational  studies  of  the  past 
30  years,  yet  no  satisfactory  theory  exists  for  turbulent  boundary  layers  with 
strong  pressure  gradients.  The  fundamental  difficulty,  to  be  sure,  stees  from 
the  general  closure  problem  for  the  time-averaged  equations.  The  recog  ltlon  of 
this  difficulty  has  led  to  the  development  of  semiemplrical  model  equations  of 
the  turbulence  and  computation  thereof,  especially  since  the  publication  of  the 
two-volume  proceedings  of  the  1968  AFOSR-IFP  Stanford  Conference  (Kline  et  al. 
1969;  Coles  6  Hirst  1969)  and  the  increasing  availability  of  larger  and  faster 
computing  resources.  This  effort  apart,  the  theoretical  approaches  resting  on 
firmer  grounds  for  addressing  pressure-gradient  boundary  layers  have  been 
twofold.  The  first  is  based  upon  similarity  analyses  and  dimensional  arguments 
(see,  e.g.,  Yaglom  1979).  The  second  is  based  upon  a  systematic  treatment  by 
the  method  of  matched  asymptotic  expansions  (see,  e.g..  Van  Dyke  1975;  Kevorkian 
ft  Cole  1981).  The  present  paper  Is  concerned  with  exploiting  this  latter 
approach,  in  the  limit  of  large  (turbulent)  Reynolds  number,  to  study  a 
turbulent  boundary  layer  under  a  strong  adverse  pressure  gradient  In  a  steady, 
two-dimensional,  incompressible  flow  past  a  smooth,  nonporous  wall. 

The  general  adverse-pressure-gradient  problem,  in  equilibrium  or  self- 
similar  form,  initially  conjectured  by  Clauser  (1954),  from  an  analogy  with  the 
zero-pressure-gradient  case,  has  been  analyzed  in  detail  by  Townsend  ( 1956a, b, 
1960,  1961a, b)  and  Rotta  (1962).  Stratford  (1959a, b)  first  studied,  both 
theoretically  and  experimentally,  boundary  layers  held  at  incipient  separation, 
the  limiting  or  strong  adverse-pressure-gradient  case  for  equilibrium  layers. 
Further  discussions  of  the  theoretical  and  experimental  results  and  their 
Implications  for  this  limiting-case  problem  are  summarized  in  Clauser  (1956), 

Mel lor  ft  Gibson  (1966),  Mellor  (1966),  Townsend  (1976),  Kader  ft  Yaglom  (1978). 
Schofield  (1981),  and  Afzal  (1983),  among  others.  Indeed,  the  results  of  this 
paper  are  compared  with  those  of  Kader  ft  Yaglom,  obtained  from  similarity  and 


dimensional  arguments,  and  with  those  of  Afzal,  obtained  from  asymptotic 
analysis,  as  well  as  with  the  original  results  of  Stratford. 

The  asymptotic  analysis  herein  for  the  "sharply  rising  pressure"  case  of 
Stratford  is  based  upon  the  Reynolds  time-averaged  equations  for  steady  two- 
dimensional  flow.  For  this  case,  as  seen  in  figure  1,  in  addition  to  the 
exterior  (or  invlscid-flow)  region,  there  is  the  Interior  (or  turbulent 
boundary-layer-flow)  region  that  consists  of  two  principal  layers.  The  first 
principal  layer  is  the  relatively  thick  nondefect  outer  layer,  characterized,  to 
leading  order  of  approximation,  by  a  convection — pressure-gradient — turbulent- 
stress  momentum  balance.  The  second  principal  layer  is  the  relatively  thin 
inner  layer,  in  which,  to  leading  order  of  approximation,  the  pressure  gradient 
and  the  turbulent  and  viscous  stresses  constitute  the  momentum  balance. 

In  Section  2,  the  steady,  two-dimensional  equations  of  continuity  and  mean 
transfer  of  x-  and  y-  momentum  and  the  complementary  boundary  conditions  for  the 
problem  are  presented  (see  Yajnik  1970;  Mellor  1972).  In  Sections  3  and  4,  the 
asymptotic  analyses  for  the  flow  in  the  exterior  invlscid  and  interior  viscous 
(or  boundary- layer)  regions  are  presented,  once  the  orderings  of  the  pertinent 
parameters  are  identified. 

For  the  boundary- layer  region,  the  postulated  characterizations  of  the  flow 
variables  in  the  outer  and  inner  layers  are  presented  in  Sections  4.1  and  4.2, 
respectively.  When  it  is  taken  that  the  streamwlse  velocity  has  a  square-root 
variation  with  respect  to  the  normal  coordinate  in  an  arbitrary  "overlap"  region 
intermediate  to  the  outer  and  inner  layers,  and  that  the  Reynolds  stress  has  a 
corresponding  linear  variation,  as  is  observed  experimentally  (see,  e.g., 
Stratford  I989a,b),  the  matching  of  the  postulated  zeroth-order  outer- layer  and 
inner-layer  solutions  for  these  flow  quantltltes  is  demonstrated  (see  Section 
4.3).  Indeed,  this  matching  provides  the  inner-edge  boundary  conditions  for  the 
higher-order  outer-layer  boundary- value  problems  and  the  outer-edge  boundary 
conditions  for  the  higher-order  inner-layer  boundary-value  problems.  In  the 
analysis  of  Section  4.1,  a  model  of  the  turbulent  length-scale  function, 
consistent  with  the  experimentally  determined  nearfleld  behavior  of  the  outer- 
layer  flow  field  and  with  the  mixing-length  hypothesis  (for  the  eddy 


dlffusivity) ,  is  developed,  as  is  a  aodel  for  the  eddy  viscosity  in  the  farfield 
of  the  outer  layer  (see,  e.g.,  Clauser  1956;  Mellor  &  Gibson  1966;  Mellor  1966). 

In  Sections  5.1  and  5.2,  respectively,  similarity  formulations  for  the 
zeroth-order  approximations  of  the  asymptotic  representations  for  the  outer  and 
inner  layers  are  presented.  For  previous  work  on  such  zeroth-order  similarity 
formulations  for  both  the  outer  and  the  inner  layers,  see,  e.g.,  Townsend 
(1976),  Mellor  6  Gibson  (1966),  Mellor  (1966),  and  Afzal  (1983). 

For  the  zeroth-order  approximation  for  the  outer  layer,  subject  to  the 
introduction  of  the  aforementioned  model  closure,  which  Incorporates  the 
experimentally  observed  nearfleld  and  farfield  behaviors  of  the  streaawlse 
velocity  and  Reynolds  stress,  the  similarity  boundary-value  problem  becomes  an 
eigenvalue  problem  for  the  determination  of  HQ,  the  (constant)  leading-order 
shape  factor,  where  this  shape  factor,  Hq.  is  related  to  the  leading-order 
streamwlse-preasure  distribution,  p^(x),  by 

HQ  -  [{l-2p*(x)}{p*'(x)}'2{-p*  (x) }-4]  -  const. 

Thus,  the  similarity  formulation  holds  for  only  a  specific  pressure 

•  • 
distribution,  pQ(x),  and/or  outer-edge  streamwise-veloclty  distribution,  uQ(x), 

i.e. . 


Pq(k)  -  2/(H0+2)]  and/or 


“o'*»  •  <*  -  ■  “oo'”  -  V'l/(,°‘a) 


This  velocity  relation  was  first  obtained  by  Townsend  (1960).  The  normal-length 
function,  hQ  (x) ,  is  (classically)  linear  in  x,  i.e.. 


hj(x)  -  [{1  -  2p* ( * ) } { Pq  (x))_1l  -  <«0  ♦  2) (x  -  xQ). 


[3] 


As  to  the  derivatives  of  this  pressure  distribution,  similarity  requires 


•2 

.  “OP  _  »“(Hq+4)/ (H.+2) 

P0  (x)  (Hq+  2)  (X  X0)  °  °  >  °’ 


*« 


«  ■  -  IBp-4l5°  u  -  .0>-2(V3>/,V2>  <  o . 

(Hq  ♦  2)  0 


With  pQ  (x)  <  0.  the  present  similarity  formulation  holds  for  a  decreasing 
adverse  pressure  gradient,  rather  than  an  increasing  one,  as  studied  by  Samuel  A 
Joubert  (1974).  Further,  p.  (x)  <  0  violates  the  moving-equilibrium  assumption 

o  •  > 

(Kader  A  Yaglom  1978;  Yaglom  1979),  wherein  it  is  taken  that  pQ  (x)  varies  only 
slowly  with  x. 


The  inner-layer  zeroth-order-approxiaation  similarity  formulation  follows 
without  the  adoption  of  a  closure  hypothesis,  in  that  only  normal  derivatives  of 
unknown  functions  appear  explicitly  in  the  zeroth-order  problem  of  Section  4.2. 
However,  this  similarity  formulation  is  based  upon  a  square-root  variation  of 
the  streaawlse  velocity  (with  respect  to  the  appropriately  scaled  normal 
coordinate)  and  a  linear  one  of  the  Reynolds  stress  at  the  outer  edge  of  this 
inner  layer.  These  variations  are  consistent  with  the  asymptotic  behavior  of  an 
inner-layer  mlxing-length-hypothesis  closure.  To  display  the  detailed  farfleld 
and  nearfleld  behaviors  of  the  flow  quantities,  here,  an  inner-layer 
dissipation-factor  closure  model  (see,  e.g.,  Van  Driest  1956;  Patankar  A 
Spalding  1970;  Szablewski  1970)  is  adopted.  From  the  present  formulation,  the 
characteristic  streamwise-veloclty  and  normal-length  functions  for  this  inner 
layer  are,  respectively. 


v*>  -  <p>»,/3  -  <*  -  x#)-«V4)/a«V8» 


*>>  -  <p*W/3-  (» -  x0,«V‘»/3'V2' 


Since  the  matching  of  the  zeroth-order  outer-  and  inner-layer  solutions  is 
established  (in  Section  4.3)  without  the  requirement  of  similarity,  the  matching 
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also  holds  upon  the  Introduction  of  similarity.  In  the  analyses  of  Townsend 
(1076),  Nellor  A  Gibson  (1966),  and  Mellor  (1966),  the  corresponding  outer-  and 
inner-layer  similarity  solutions  are  patched  (at  a  point  of  assumed  common 
validity),  rather  than  matched.  In  the  analysis  of  Afzal  (1983),  the  zeroth- 
order  inner-layer  similarity  boundary-value  problem  is  essentially  correct; 
however,  because  the  matching  requirement  for  the  streamwlse  velocity  is  not 
correctly  formulated,  the  given  complementary  outer-layer  similarity  boundary- 
value  problem  is  subject  to  incorrect  inner-edge  boundary  conditions.  Afzal 
matches  the  (appropriately  scaled)  normal  derivative  of  the  inner-edge  limit  of 
the  outer-layer  streamwise-velocity  solution  and  that  of  the  outer-edge  limit  of 
the  inner-layer  streamwise-velocity  solution.  The  correct  procedure  (see.  e.g., 
Kevorkian  A  Cole  1981)  requires  the  direct  matching  of  the  inner-edge  limit  of 
the  outer-layer  streamwise-velocity  solution  and  the  outer-edge  limit  of  the 
inner-layer  streamwise-velocity  solution  in  an  arbitrary  intermediate  layer. 

The  same  procedure  is  required  for  the  matching  of  the  corresponding  Reynolds- 
stress  solutions.  In  the  present  case,  the  introduction  of  Stratford's 
empirical  results  and/or  the  specification  of  the  mixlng-length-hypothesls 
closure  provides  the  limiting  behaviors  of  the  outer-  and  inner-layer  solutions, 
and  the  zeroth-order  matching  is  accomplished.  In  turn,  it  is  possible  to 
determine  the  appropriate  higher-order  terms  in  the  outer-  and  inner-layer 
asymptotic  expansions.  These  expansions  are  Introduced  in  Section  4;  their 
correctness  is  confirmed  in  Section  6. 

From  the  analyses  of  Sections  4.1  and  4.2,  the  equations  of  notion  for  the 
higher-order  approximations  for  both  the  outer  and  the  inner  layers  of  the 
interior  viscous  region  are  obtained.  From  matching  with  the  exterlor-inviscld- 
region  solutions,  the  boundary  conditions  at  the  outer  edge  of  the  outer  layer 
are  obtained;  the  boundary  conditions  at  the  inner  edge  of  the  inner  layer 
(i.e.,  the  wall)  follow  directly  from  the  original  boundary-value  problem;  and 
the  hlgher-order-approxiaation  boundary  conditions  at  the  inner  edge  of  the 
outer  layer  follow  from  matching  with  the  outer-edge  solutions  of  the  zeroth- 
order  approximation  for  the  inner  layer,  just  as  the  hlgher-order-approxiaation 
boundary  conditions  at  the  outer  edge  of  the  inner  layer  follow  from  matching 
with  the  inner-edge  solutions  of  the  zeroth-order  approximation  for  the  outer 
layer.  With  this  information  and  that  obtained  from  the  zeroth-order  outer-  and 


inner-layer  similarity  boundary- value  probleas,  it  is  possible  to  develop  a 
sequence  of  higher-order  similarity  boundary-value  problems  for  both  the  outer 
and  the  inner  layers.  The  details  of  this  development  are  presented  in  Section 
6.  It  is  shown  that  the  outer-layer  solutions  for  the  streamwlse-velocity  and 
Reynolds-stress  functions  are  of  the  following  forms: 


a_  *  ,  a./2  *1l  ,  .  .  ,3/2  "a  . 

*  dr*  *  rfn  *  dri  •  rln 


r  ,  1/2 ,  .  ,  .3/2 , 

“5F  -  *0  +  A,  *  A.  ^ 


+  .  .  .  . 


Here,  u  -  u(x,y)  and  r  -r(x.y)  are  the  (nondlmensional)  streamwlse  velocity  and 

•  • 

turbulent  shear  stress  in  the  outer  layer,  respectively;  u  ' »  u  (x)  - 

*  1/2  00  * 

(1  -  2p  q(x)]  is  the  zeroth-order  outer-edge  streamwlse  velocity;  rj  -  y/h  is 


the  outer-layer  similarity  coordinate,  where  y  is  the  appropriately  scaled 

•  *  *  *i  _i 

outer-layer  normal  coordinate  and  hQ  »  hQ(x)  -  [{1  -  2pQ(x)}{p0  (x)}  ]  is  the 


normal-length  function  for  this  layer;  fk  -  fk(v)  and  ^  tj )  are  kth- 
approxlmatlon  similarity  streamf unction  and  Reynolds  stress,  respectively;  and 


A,  -  A .<«:<)  -  «t(i  -  2p’(«))'‘(p*'(x))a/3] 


.  3K.„ .  2)  »•  (« - .  ,(vi»/<v*>r,/*. 


-  -  ~  ~  -9/7 

where  6  «  (6/ c)  «  ( ujc/v )  is  the  nondlmensional  thickness  of  this  outer 


layer.  It  is  also  shown  that  the  corresponding  inner-layer  solutions  for  the 
streamwlse-velocity  and  Reynolds-stress  functions  are  of  the  following  forms: 


y  *  3?  • 


*  *0  +  A*  *2  +  • • * 


VV.V.V.V.VA 


Here,  u  -  u(x,r)  and  r  -  r(x,r)  are  the  (nondiaensional)  streaawise  velocity  and 

•  • 

turbulent  shear  stress  In  the  inner  layer,  respectively;  u  .  -  u  (x;d)  - 

1/2  ••  1/3  p0  p0 

6  (p  ( x ) >  is  the  (so-called)  pressure  velocity,  the  characteristic  zeroth- 

0  * 
order  streaawise  velocity  of  the  inner  layer;  C  »  r/hp0  is  the  inner-layer 

siallarlty  coordinate,  where  r  -  y/6  is  the  inner-layer  normal  coordinate  and 
•  *  *'  -1/3 

h  _  •  h  .  (x)  -  {p.  (x)}  is  the  normal-length  function  for  this  layer;  and 
po  po  0 

Kk  "  8^(0  And  ”  «*^(  C)  are  the  kth-approxlaation  siallarlty  streaafunctlon 
and  Reynolds  stress,  respectively. 

Two  important  results  are  obtained  froa  the  siallarlty  formulations  of 
Section  6:  (1)  the  siallarlty  coordinates,  17  and  C,  respectively,  are  related, 

through  the  expansion  paraaeter,  a,,  by 


-1 


jj  -  and/or  C  -  a,  ij 


•  * 

and  (2),  the  characteristic  streaawise  velocltites,  uQ  and  upQ,  are  also 
related,  through  the  expansion  paraaeter,  a,,  by 


*  1/2  *  *  -1/2  * 

upo  "  A*  uo  and/or  uo  -  A*  u 


pO 


In  Section  7.  the  solutions  for  the  streaawise  velocity  and  the  Reynolds 
stress  are  presented  for  the  distinguished  similarity  intermediate  layer,  the 
existence  of  which  is  suggested  by  the  inner-edge  behaviors  of  the  outer-layer 
expansions  and/or  the  outer-edge  behaviors  of  the  inner-layer  expansions,  found 
in  Section  6.  For  this  interaediate  layer,  the  siallarlty  normal  coordinate  is 


.  ..1/2  -1/2 

( 17C)  -  A.  T) 


'i/2< 


For  a  discussion  of  three-layer,  or  intermediate-layer,  theory,  within  the 
fraaework  of  (classical)  two-layer  theory  for  turbulent  pipe/channel  and 
favorable-pressure-gradient  boundary- layer  flows,  see,  e.g.,  Afzal  (1982);  Afzal 


and  the  appropriate  speed  function  is 


*  *  •  1/2 

-  (urtu^)1/2 


u0  “  A* 


In  turn,  the  lnteraediate-layer  solutions  for  the  streaawlse-velocity  and  the 
Reynolds-stress  functions  are 


_JL_  oi  ra*  „1/2i  *  A1/4ra+  i 
•  -  W01x  1  •  [^in] 


10J 


‘10 


.  a1/2 r-+  v3/2  +  -1/2.  .3/4,.+  .  +  -1,  . 


21' 


12*  32 


k~  «  *  o  *  ^/4t*;a«3/2  <■  ->;,x'i/2i * ... . 


u 


10 


where  0*^’  and  2*^  are  consts.,  whose  values  depend  upon  the  outer- 


and  inner-layer  closures  employed.  Consistent  with  Stratford's  data  and  the 
nixing-length-hypothesis  closure. 


0qi  ”  2 /k  (*  4.9  for  k  «  0.41). 


The  representations  of  Kader  A  Yagloa  (1978),  Afzal  (1983),  and  Stratford 
(1959b)  for  the  experimental  data  on  the  streaawise  velocity  in  the  lnteraediate 
layer  are  re-cast  in  the  forw  of  the  aforaentioned  interaediate-layer 
streaawise-velocity  solutions.  The  results  Indicate  that  predictions  based  on 
the  three-layer  theory  are  better  than  those  based  on  the  (classical)  two-layer 
theory . 
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2.  EQUATIONS  OP  MBAM  NOTION 

! 

Consider  the  steady  two-dimensional  fully  developed  turbulent  flow  of  an 
incompressible  fluid  of  constant  density  and  (kinematic)  viscosity  (p,  v  - 
consts.)  past  a  smooth  flat  surface  (cf.  Yajnik  1970;  Mellor  1972).  Let 
X  -  c  X  and  Y  =  c  Y  represent  the  coordinates  tangential  and  normal  to  the 
surface,  respectively,  with  c  a  characteristic  body  length.  The  mean  velocity 
components  in  the  X-  and  Y-directlons ,  respectively,  and  the  mean  pressure  are 
«  uJJ,  V  ■  u^V,  P  *  p^  +  pu2  P,  with  uw  and  pw  being  the  velocity  and  pressure 

in  the  undisturbed  region  far  from  the  surface.  The  laminar  stress  tensor 

~  — 2  ~  ~  ~  ~ 

components  are  Sjj  ■  pu^u^/u^Sj  j  (i  or  J  -  X,  Y  and/or  X,  Y),  with 
S..  -  9U./3X.  and  u  -  v/c,  a  characteristic  "laminar  velocity";  the  turbulent 

*  J  "  m  - 2  w  M  2  m.  mm 

stress  tensor  components  are  T. .  *  pu  (u./u  )  T. .  (i  or  j  -  X,  Y  and/or  X,  Y) , 

-*  lj  •  t  *  lj 

with  ut  a  characteristic  "turbulent  velocity.”  In  the  analysis  that  follows,  it 
is  taken  that  the  turbulent  and  laminar  velocity  parameters,  e  =*  (u./u  ) 

-  ~  -  -i 

and  p  =  (ut/uw)  -  (v/ujc)  -  R  ,  respectively,  go  to  zero. 


In  the  domain  X  >  XQ,  0  <  Y  <  the  (nondlmenslonal)  equations  of  mean 
motion  for  this  flow  are: 


au  av  .  „  aw  ,,  aw 

ax  *  ay  *  0:  u  “  3Y’  v  "  "  ax  : 


(2.1) 


(uau  +  v^U)  +  3P  .  e2  (!!m  +  !lxv 
'  <ix  oY'  ax  'ax  av 


,a2u  a2u. 

ax  ayz 


(2.2a) 


(02Y  ,  yM,  ,  JE  .  e2  (iLa  ,  fLoc 

'  cix  ^y'  ay  'ax  ay 


,3^  aV 


)/  W  »  V  »  i 

+  «  2  *  ~2] 
ax  aYc 


(2.2b) 


The  surface  boundary  conditions  are: 


U,  V  -  o,  Ttj  -  0  as  Y  -  0  (X  >  XQ); 


(2.3a) 


the  freestreaa  boundary  conditions  are: 


a  i 


The  upatreaa  boundary  conditions  are  taken  to  be 


X,  Y  fixed  as  t,  n,  and  d(e.jz)  ■*  0; 


(2.4) 


those  in  the  two  principal  layers  of  the  interior  region,  respectively,  are 
studied  for 


$ 


5; 


U.  V,  P,  Tjj  -  fncs(Y)  as  X  -  XQ  (0  <  Y  <  «). 


(2.3c) 


In  the  present  analysis,  the  existence  of  an  exterior  (or  inviscld-flow) 
region  and  an  interior  (or  turbulent  boundary-layer-flow)  region  is  stipulated. 
Further,  it  is  taken  that  this  boundary-layer  region  consists  of  two  principal 
layers.  By  aeans  of  liait-process-expansion  techniques,  the  solutions  for  the 
flow  quantities  in  the  exterior  region  are  studied  for 


I 


x,  y  fixed,  with  x  =  X,  y  =  Y/3.  as  e,  //,  and  8(e,/i)  ■*  0; 


x,  r  fixed,  with  x  =  X,  r  *  Y/8  ,  as  e,  n,  and  6(e,n)  ■*  0. 


Here,  8(&,n)  is  the  effective-thickness  paraaeter  of  the  interior  region,  as 

2 

well  as  that  of  the  outer  layer  of  this  region;  8  (e,jt)  is  the  effective¬ 
thickness  paraaeter  of  the  inner  layer  of  this  interior  region. 

In  this  study,  it  is  deteralned  that  e  and  n  are  not  independent 
paraaeters;  rather,  it  is  deteralned  that  e  =  e(p),  and,  in  turn,  8{e,n)  *  8(n) . 
The  following  analysis  requires  that 


(2.5a) 

-  :? 
0.  -V 

«v  v 

(2.5b) 

'  «£ 

e(d)  -  81/Z,  and  n(8)  -  <J7/2, 


(2.6a) 


that  is. 


e(^)  -  fi1/7 ,  and  8(n)  -  Z/1 . 


(2.6b) 
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3.  THI  EXTERIOR  REGION 


For  X,  Y  fixed,  with  e,  ft  -  0,  the  exterior-region  expansions  for  the  mean- 
velocity  components,  U,  V,  and  the  pressure,  P,  are  in  integral  powers  of  the 
boundary-layer-thickness  parameter,  8  -  d(e,ft),  namely: 

G(X,  Y;  «.  ft)  -  G(X,  Y;  8)  £  GQ(X.  Y)  ♦  dG^X.  Y)  *  d2G2(X,  Y)  +  ...  ,  (3.1) 

where  G  -  U.  V,  P.  With  the  stipulation  that  there  is  no  exterior-region 
turbulence  and  that  the  interior-region  turbulence  decays  exponentially  as  the 
exterior  region  is  approached,  it  is  taken  that  the  exterior-region  expansions 
for  the  turbulent-stress  components  are  in  transcendental ly  small  powers  of  8, 
such  that 

TjjfX.  Y;  e,  fx)  -  ^.(X.  Y;  8)  ~  0.  (3.2) 

Substitution  of  the  expansions  of  (3.1)  and  (3.2)  into  (2.1)  and  (2.2) 
yields  equations  for  the  successive  exterior-region  approximations.  The 
equations  for  the  zeroth-,  first-,  and  second-order  approximations, 
respectively,  are 


3U0  avo 

a¥o 

a*0 

dX  +  3Y 

-  0:  U 

"  3Y 

.  v  -  -  — -  ’ 

■  o  ax 

(3.3a) 

3U0 

auo 

apo 

(uo  ax  + 

v  — -) 

f  ax 

o. 

(3.3b) 

avo 

dVo 

apo 

(uo  ax"  + 

v  — -) 

0  3Y  ' 

f  W  m 

0; 

(3.3c) 

OHi  avi 
ax  +  ay 

afi 

avl 

-  0:  Uj 

"  3Y-’ 

v - - 

i  ax  • 

(3.4a) 

au 

au 

auo 

au  ap 

(uoix  + 

V0  3Y-  + 

aTui 

+  3Y  V  +  ax"  "  °’ 

(3.4b) 

avi 

avl 

avo 

avo  3P1 

(uo  ax  + 

V0  3Y  + 

iirui 

+  W  V  +  W  ”  0: 

(3.4c) 
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aU2  3U0 
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3Y  +  ax 

*  3Y 

ax 

av„ 

av„ 

dVn 

avrt 

ap„ 

i 

_ 2 

0 

U2 

* _ 2. 

V 

*  — 1 

'o 

ax 

*  vo 

3Y  ax 

*  3Y 

3Y 

3U 


au. 


-  (u 


1  ♦  V,  TT^)  ,  (3.5b) 


1  ax  1  3Y 


avl  avl 

(ui  ax  *  vi  aY 


(3.3c) 


Por  these  equations,  the  freestreaa  boundary  conditions  are 


UQ  -  1.  VQ,  P0  -  0,  Ur  Vj,  Px  -  0.  U2.  V2,  P2  -  0  as  Y  -*  •;  (3.6a) 


the  (only  relevant)  surface  boundary  condition  is 


V0  -  0  as  Y  -  0. 


(3.6b) 


In  turn,  it  is  possible  to  construct  first  integrals  of  (3.3)  -  (3.0).  Here, 
the  first  integrals  for  the  total  pressure  and  the  vorticlty  are 


- 


tpo*i(02*vS»  -  w- 


av0  au0 

(aT’ar)  “-Bo<V: 


(3.7a) 

(3.7b) 


"l  -  [P1  +  (Vi  +  V0Vi)]  “  IflB0(V  *  W1* 


(3.8a) 


avi  a«i 

Ji  Wo>  *  Bi(V1: 


(3.8b) 


n,  -  tP„  ♦  (u„u.  ♦  vftv,)  *  i(u?  ♦  v?) ] 
2  2  02  0  2  21  1 


•  ‘wv  *  w1  *  vivi'v  * 


(3.9a) 


3V2  SO, 


°2  •  (3X  '  if' 1  '  -  'VS'V  *  Bi<V>  -  V«Va"<V  *  Bl"(Vl-  13  <’b> 


With  the  freestreaa  taken  to  be  independent  of  d(e,  ft)  [cf.  (3.6a)],  it  is 
consistent  to  take  the  higher-order  Bernoulli  functions  to  be  zero,  that  is, 
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(3.10a) 


B1<V'  W' 


Further,  if  the  exterior-region  flow  Is  lrrotational  (i.e.,  -  o,  k  -  0,  1, 


2,  . ..),  It  follows  that 


VV  -i'  "dBo'V-  B5<V-  •  -«• 


(3.10b) 


Subject  to  the  constraints  of  (3.10),  the  equations  for  the  exterior  region 


can  be  written  as 


av.  au  a2f.  a2? 

(^T  -  w  >  *  -  l  9  +  9  )  “  0  (k  -  0,  1.  2. 

ax  3Y  ax2  av2 


(3.11a) 


a*  2  a*  2 

i  t!  -  {(^r)  *  (“)  >1. 


4!!a!Ii  !!o  f!i. 

'ax  ax  +  3Y  3Y  *  * 


dW0  a*2  3f0  a,2  1  afl  2  9V1  2 

-  KaraT‘irar>  ns1*  +  >3*  • 


(3.1lb) 


It  is  taken  that  the  solutions  to  the  differential  equations,  either  in  the 
forws  given  by  (3.3)  -  (3.3)  or  in  those  given  by  (3.11),  are  analytic  as  Y  -*  0. 
Thus,  near  the  surface  (Y  •*  0),  the  exterlor-reglon-flow  quantities,  6  -  U,  V, 


P,  can  be  written  as 


.  30  •  3*0  . 

t°0  *  'av2'  y  *  y  *  -  l 


**1  e  2  . 

♦  «  [Oj  ♦  (gy1)  Y  ♦  ...)  ♦  «  (Qg  *  ...]  ♦  . 


loo  *  eo.v  *  &2y*  *  ■>  *4t°!  *°!iy* 


*  a2  (0*  .  , 


(3.12) 
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with  G^,  Gk|  -  fncs(X).  The  surface  speed  UQ  arises  from  the  basic  inviscld 
flow  [of  (3.3)];  the  surface  boundary  condition  [of  (3.6b)]  for  this  basic 
inviscld  flow  gives  VQ  -  0.  Based  upon  these  equations  of  aotlon.  it  is 
determined  that 


11  ’  dX  ’ 


dU0  • 
dF'  V02  "  °* 


*  1  *2  * 

P0  -  2  (1-U0  >•  P01 


dX . 


«  • 

"  U0  Ul* 


*\r %^<v*  %*  •,'  s'*,* 


(3.13a) 


(3.13b)  Q 


(3.13c) 
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4.  THE  TURBULEHT  BOURDARY  LAYER 


For  the  turbulent  boundary  layer,  the  coordinates  are 


x  -  X.  y  -  Y/4,  with  4-0; 


(4.1) 


the  flow  quantities  are 


o  -  u.  v  -  av.  p  -  p.  tu  -  t4J  . 


(4.2) 


Here.  4  -  4(e,*i);  G  -  G(X,Y;  e.n)  -  G(X,Y;  8).  with  G  -  U,  V,  P.  T^; 


S  *  g(x,y;  4),  with  g  >  u,  v,  p,  t^y  Introduction  of  these  boundary-layer 
variables  into  (2.1)  and  (2.2)  yields 


3u  dv  ^ 

to  *  iP  “  0:  u  "  ip  *  v 


M 


_  M  . 

ax  ' 


(4.3) 


(u  lx  +  v  lx,  +  la  „  +  s  fl**>  +  (U-u&l  *  a2  1?x» 

(u  ax  +  v  ay'  +  ax  ^'ay  +  3  ax  »  *  2  +  *  TT  * 

o  oy  ox 


(4.4a) 


2  at. 


at 


.2/..  lx  av  aj>  j|/AL  fwvyy  .  .  wtvx.  .  ,2.n  .  ,a2v  .  .2  a2v 

*  (u  ax  *  v  ay)  *  ay  a{8>  ay  +  3  to  >  *  a  (^(^2  *  3 


).  (4.4b) 


The  surface  boundary  conditions  are 


u.  v  -  0.  t4J  -  0  as  y  -  0. 


(4.5a) 


The  Interior-region — exterior-region  Batching  conditions  are 


u  -  U.  v  -  V/4,  p  -  P.  t4J-  TlJ  -  0  a.  y  -  Y  -  0. 


(4.5b) 


4 . 1  The  Outer  Layer 


The  analysis  of  the  boundary  layer  is  Initiated  with  the  consideration  of 
an  outer  layer,  defined  by  x,  y  fixed,  with  s,  u,  and  8(t,n)  -  0.  For  this 
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outer  layer,  it  la  taken  that,  to  leading  order  of  approxiaatlon,  the  outer- 
layer  x-aoemntua  equation,  (4.4a),  is  characterized  by  a  convection— pressure- 
fradient—turbulent-strees  balance  (with  the  viscous  stress  negligible). 
Specifically,  it  is  taken  that  the  outer-layer  thickness  paraaeter  is  equal  to 
the  square  of  the  turbulent -velocity  paraaeter,  that  is, 

8i9.lt)  -  8(9(11))  -  92(v)  *  0:  9(8)  -  81/2  -  0.  (4.6) 

7/2  2 

with  9(8)  -  8  -  0  (  to  be  Justified  in  Section  4.2),  such  that  (9/8  )  - 

3/2 

8  ■*  0.  Por  these  orderings  of  the  paraaeters,  (4.4a)  and  4.4b),  eaploylng 

the  thickness  paraaeter,  8,  throughout,  take  the  foras 


(4.4a)' 


ay  ay 
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*2<«  S  * v  &  ♦  5  -  *  ♦  *in  CM  ♦  *2  M)  • 


St, 

ax 


(4.4b)’ 


ay* 


ax 


In  turn,  froa  an  exaaination  of  (4.4)’,  it  is  taken  that  the  outer-layer 
expansions  for  the  dependent  variables  are  of  the  fora 


g(x.y :«)  3! 


s0(*.y) 


a1/2gj(x.y) 


ag2(*.y) 


a3/2g3(*.y) 


(4.7) 


where  g  -  u,  v,  p.  r 
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Thus,  the  leading-order  approx 1 nations  for  the  outer  layer  are 
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(4.8a-c) 


(4.Sa-c) 
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'  o  ax  o  ay  ax  2  ay  2'  ax  ay 
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,  au3  fUa  3u0  3u0  %  ap3  !!xy3 
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In  what  follows,  the  notation  r  -  r  -  t,  and  r  ,  •  r  ,  -  r,  , 

xy  yx  xyk  yxk  k 

k  -  0,  1,  2,  3 . la  employed. 


In  this  paper,  for  specificity,  the  slxlng-length-theory  closure  Is  adopted 
for  the  nearfleld  behavior  of  the  outer  layer.  That  is.  the  turbulent  length- 
scale  function,  I,  is  taken  to  be 


I  -  r  <"ay )  1  ■*  *Y  y  -  0.  "1th  x  -  von  Karwan  const. 


(4.12) 


In  terws  of  the  outer-layer  expansions,  this  nearfleld  behavior  (y  -»  0)  can  be 

expressed  as 
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Thus,  the  zeroth-order  approximation  for  the  outer  layer  is 
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(4.14a) 

(4.14b,c) 


The  outer-edge  (y  -*  •)  boundary  conditions  are 
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where  uQ  and  pQ  are  the  externally  applied  streanwise  velocity  and  pressure, 
respectively.  Based  on  an  external  flow  that  is  lrrotatlonal  (cf.  Section  3), 


*  1  *2 

it  follows  that,  to  leading  order  of  approximation,  p  -  -(l-u  ),  such  that 

*  *  •  u  z  o 


(dpQ/dx)  «  -  uQ(duQ/dx)  >  0  tor  the  adverse-pressure-gradient  case  under 


consideration  here.  The  (only  relevant)  surface  (y  *  0)  boundary  condition  is 
vQ  ”  0.  (4.15) 


Directly,  it  is  seen,  from  (4.14c),  that  pQ  -  pQ  -  fnc(x).  In  turn,  (4.14b) 


becomes 
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'  “n  ts  a.,  f  a*.  a*. 


(4.14b) 


'”0  dx  0  ay  '  dx  dy 

Based  upon  the  preceding,  the  following  integral  relations  are  developed: 
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for  rQ  ■*  v+  as  y  ■*  0,  with 
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"J-)  dy. 


(4. 18a,b) 


the  zeroth-order  defect  and  aoaentua  thicknesses,  respectively. 


Based  on  the  results  of  Stratford  (1959a),  it  Is  anticipated  that,  as 
1/2 

y  ■*  0,  uQ«  y  -*  0  and  y  •*  0.  Specifically,  the  solutions  of  (4.14a)  and 
(4.14b)',  as  y  -*  0,  taking  into  account  the  nearfield  closure  hypothesis  of 
(4.8d),  are 
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These  inner-edge  (y.jj  -  0)  results  nay  be  expressed  as 
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Further,  the  outer-edge  (y.ri  *  «•)  results  say  now  be  expressed  as 
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4.2  The  Inner  Layer 

The  results  of  (4.19)  and  (4.20)  Indicate  that  the  outer-layer  flow 

variables  u  and  r  ■*  0  and  (3u  /3y)  -  •  in  the  Unit  of  y  -*  0,  3  -*  0;  however, 

9  9  1/2  0 

It  is  seen  that  uQ/d  and  tQ/3  are  fixed  for  r  -  y/3  fixed.  Based  on  these 
results,  it  is  appropriate  to  introduce  a  layer  interior  to  the  outer  layer.  It 
is  in  this  inner  layer  that  the  surface  boundary  conditions  are  satisfied 
unlforaly.  Based  on  the  aforementioned  results,  this  inner  layer  is  defined  by 


(21] 


(4.24) 


x  -  X,  r  -  y /8  -  Y /8  ,  with  8  ■*  0, 


with  the  flow  quantities  given  by 


U  =  u  =  81/Zm,  V  =  8v  =  d5/2n,  P  -  p  -  q.  =  St^  .  (4.25) 


Here,  G  =■  G(X,Y;8),  with  G  *  U,  V,  P,  T^;  g  =  g(x,y;3),  with  g  -  u,  v,  p, 


f  -  f(x,r;6),  with  f  =  »,  n,  q,  tjj.  In  terws  of  the  inner-layer  variables,  the 


equations  of  Motion  are 
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In  order  that,  to  leading  order  of  appr ox i nation,  the  viscous-stress 
contribution  be  of  the  sane  order  of  Magnitude  as  the  turbulent-stress  and 
pressure-gradient  contributions  in  the  x-aocsntun  equation,  (4.27a),  it  is 
necessary  that 
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n(8)  »  8  -*  0  as  8  •*  0, 


(4.28) 


1/2  -1  ~  -  ~ 
as  well  as  e(3)  >3  -  0  [see  (4.6)].  Thus,  since  n  **  R  ,  where  R  =  (u  c/v) 


is  the  reference  Reynolds  nunber,  it  is  deternined  that  the  present  analysis 
holds  for 


3(R)  -  R~2/7  -  0,  e(R)  -  R_1/7  -  0  as  R  -  .. 


(4.29a,b) 


Por  this  ordering  of  the  paraaeters,  (4.27a,b)  can  be  written  as 
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Proa  an  examination  of  (4.26)  and  (4.27)',  the  inner-layer  expansions  for  the 
dependent  variables  are 
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f (x,r ;d)  £  fQ(x,r)  ♦  a1/2f1(x,r)  +  3f2(x,r)  +  d3/2f3(x,r)  ♦  ...  ,  (4.30) 


where  f  »  a,  n,  q,  t^j. 
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The  resulting  equations  for  the  leading-order  approximation  for  this  inner 
layer  are 
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(4.34a-c) 
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In  what  follows,  the  notation  t  -  t  -  t  and  t  ,  »  t  ,  -  t,  , 

xr  rx  xrk  rxk  k 

k  -  0,  1,  2,  3 . Is  eaployed. 

Thus,  the  equations  for  the  zeroth-order  approximation  are 
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0)  boundary  conditions 
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In  this  approximation,  from  (4.36b),  the  pressure  is  found  to  be 

Q0  *  q0  “  fnc(x)  •  (4.38) 

Directly,  the  first  integral  of  the  x-momentum  equation,  (4.36a),  is  determined 
to  be 


0  +  dr  dx  r 
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(4.39) 


For  the  inner  layer,  the  appropriate  turbulent  length-scale  function  is  A, 
defined  by 
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At  the  outer  edge  (r  -*  »)  of  the  inner  layer,  the  mlxing-length-theory  closure 
of  Aq  *  kt  -  •»  holds.  This  behavior  is  consistent  with  the  previously  adopted 
inner-edge  behavior  of  the  outer-layer  length-scale  function  [i.e.,  t  ~  ry  -*  0 
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as  y  -*  0].  For  AQ/*r  ■*  1  (exponentially)  m  r  «  •,  from  (4.39)  and  (4.40),  It 
follows  that,  with  q*  «  q®  and  t*  =  (3m0/ar)°, 
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Integration  produces 
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where  a0.  ...  -  fncs(x)  (to  be  determined).  In  turn, 
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These  outer-edge  results  may  also  be  expressed  as 
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4.3  Matching  of  Zeroth-Order  Solutions 


A  cowparlson  of  (4.19),  (4.20),  and  (4.42),  (4.43)  Indicates  that  the 

zeroth-order  outer-layer  and  inner-layer  solutions  for  the  streaawlse  velocity 
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Purther  consideration  of  the  watching  of  the  outer-layer  and  inner-layer 
solutions  yields 
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5.  ZEROTH -ORDRRAPPROXINAT ION  SIMILARITY  FORMULATIONS 


5.1  The  Zeroth-Order  Approximation  for  the  Outer  Layer 


The  results  of  (4.20)  -  (4.23)  suggest  that  a  similarity  formulation  for 

this  zeroth-order  approximation  exists.  For  this  similarity  formulation,  with 

•  *  1/2  *  *  *  -1 

uQ  *  [(l-2pg)  ] ,  hjj  *  [ ( 1-2Pq) (dpg/dx)  ],  the  independent  variables  are 
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and  the  dependent  variables  are 
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Proa  the  nearfleld  mixing-length  hypothesis,  similarity  requires 
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This  requirement  with  respect  to  A  indicates  that  similarity  holds  only  for  a 

* 

specific  class  of  applied  pressure  distributions,  pQ(x).  More  will  be  said 
about  this  later. 


Introduction  of  these  similarity  variables  into  the  equations  of  notion  for 


this  zeroth  approximation  produces 
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The  farfield  boundary  conditions  are 
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the  nearfleld  boundary  conditions  are 
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Integration  of  (5.4)  over  the  doaaln  of  tj  ,  subject  to  (5.5) 
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Thus,  for  this  slallarlty  foraulatlon,  the  shape  factor,  H., 
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applied  pressure  distribution,  pQ,  by 
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The  solution  of  this  equation  for  pQ  is 
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where  B  and  C  are  constants  of  Integration.  Quantities  associated  with  this 
solution  are 
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With  the  appropriate  choice  of  the  constants  B  and  C,  the  results  of  (5.10)  can 
be  written  as 
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In  terns  of  HQ.  rather  than  A,  the  sinilarity  boundary-value  problen  for 
the  zeroth-order  approxinatlon  for  the  outer  layer  is 
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5.2  The  Zeroth-Order  Approximation  for  the  Inner  Layer 

A  slallarlty  formulation  for  thia  approximation  exists.  Based  on  (4.44) 
and  (4.45),  for  this  slallarlty  foraulatlon.  the  Independent  variables  are 
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Here, 
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With  respect  to  the  similarity  length-scale  function,  0Q  -  /Q 

1/2  22-1  U  O 

•*  (d  g0/dC  )  .  the  following  asymptotic  behaviors  hold: 


1/2  d  20  * 
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#0-*i/2  <^T>  1  ~  'S.  *#<1  +  ...|  0  as  {  »  0, 

At 


with  s,  #0>,  ...  -  consts. 


These  behaviors  suggest  the  following  approximate  closure: 
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For  this  closure,  0.  *  *C(1  +  ...)  as  C  ■*  •»,  and  0fl  ~  0~  ?  (1  +  ...)  as  C  -*  0, 
with  s-2,  -  (*•  /Z).  Por  this  closure  hypothesis,  (5.16)  takes 
the  fora 


2  2  2  **  *0  ^  ^  ®o  + 

*CE  (— f)  +  (— r2)  -  (•  +0-0. 

dCZ  d?z 


(5.21) 


Thus, 
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For  C  -  0,  the  asymptotic  behaviors  of  (dg^/dO  and  »Q  are  now  determined  to  be 
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For  C  -*  «*.  the  asymptotic  behaviors  of  (dgQ/dO  and  wQ  in  this  limit,  as  given 
in  (5. 17b) ,  hold. 


Note  that,  for  the  closure  of  (5.20), 
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6.  HIGHER -ORDER -APPROXIMATIONS  SIMILARITY  FORMULATIONS 
6.1  Higher-Order  Approximations  for  the  Outer  Layer 


The  equations  of  aotion  for  gjfx.y),  with  gj  *  Uj ,  v Pj ,  have  been 


given  previously  [see  (4.9)].  Froa  consideration  of  the  aatching  with  the 
exterior  layer, 


V  pr  'iji  *  0  *•  y  -  - 


(6.1) 


Thus,  Pj  *  Pj  ”  fnc(x)  *  0.  From  aatching  of  the  outer  layer  with  the  inner 


layer  [see  (4.48)]. 


♦  +  dP0  1/3 

Uj  -  a0  -  r  (^“)  .  -  0  as  y  -  0. 


(6.2) 


+  * 

Recall  that  r  =  const,  and  (dpQ/dx)  -  fnc(x;HQ)  are  introduced  in  the  zeroth- 


order  outer-  and  inner-layer  similarity  formulations. 


For  the  first-order  outer-layer  similarity  formulation,  the  independent 
variables  are  transformed  from  (x.y)  to  (x.jj),  where  7j(x,y)  is  defined  in  (5.1). 
The  first-order  dependent  variables  are 


•  dp  -2/3 

Pjtx.y)  -  ^(X.n)  -  [(1-2P0Hdx  )  1  fj(i|)5 
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Introduction  of  (6.3)  and  (5.2)  Into  (4.9),  subject  to  (6.1)  and  (6.2),  produces 
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rf2fl 

2  d2f0 

dr, 

-  -  t  (H. 
3  0 
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(6.5a) 
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2  2 
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dr, 

-  T  .  #1  -  *«v  ~2  ~2  "  ° 

dr)  dr) 

as  i)  ■*  0. 

(6.5b) 

Proa  (6.4)  and  (6.5),  the  inner-edge  (r,  •*  0)  behaviors  of  the  velocity  and 
stress  functions  are 


dr,  fi10  *  ^12*  *  -  '*^10* 


with  B+  ■  b+  »  r+ ,  6+ 
p10  00  '  p12 
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Vhe  equations  of  notion  for  g9(x,y), 

*•  • 

presented  in  (4.10).  Since  p2  -*  pg  -  P 
integration  of  (4.10c)  yields 


with  gg  -  u2>  v2,  p2>  tjjg,  are 
«  fnc(x),  and  *  0  as  y  •*  •, 


P2  "  p2  *  CyyO  * 


In  turn,  (4.10b)  becoaes 


du 
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(u/»  *  v/i  HT~  +  +  ITT  V,)  ♦  —  - 
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0  dy  dx  2  dy 


dx  dy 
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Matching  with  the  exterior  layer  [see  (3.13)]  produces 
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(6.9) 


Similarly,  Batching  with  the  inner  layer  [see  (4.48)  and  (5.18)]  produces 
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(6.10) 


For  the  second-order  outer-layer  siBilarity  formulation,  the  independent 
variables  are  transformed  from  (x,y)  to  (x.jj),  and  the  dependent  variables  are 


*2(x.y) 


*  1/2  dp0  ~^3 
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(6.11b) 
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Note  that,  with  p2(x,y)  *  p 2(x,r))  -  [(dpQ/dx)  ]*2(r,)  -  -  W2(dpQ/dx) 

*  2/a 

ryyO(X,y)  ■  ryyO(X,,,)  “  HdpQ/dx)  /JH*2(r,)  ♦  Wg]  *  0.  and 
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The  following  boundary-value  problem  results: 
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The  inner-edge  (r)  -*  0)  behaviors  for  the  velocity  and  stress  functions  are 
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In  deterainlng  these  behaviors,  it  has  been  taken  that  ^  ~  $02»}  «•  ...  as  >j  -  0. 


The  equations  of  aotion  for  g3(x,y),  with  gg  =■  u3>  v^,  p3>  r^g,  are 
presented  in  (4.11).  Proa  consideration  of  the  watching  with  the  exterior 
layer , 


V  P3'  rij3  ■*  0  a8  y  - 


Based  on  (6.15),  integration  of  (4.11c)  yields 


p3  “  ryyl  • 


In  turn,  (4.11b)  becoaes 
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Matching  with  the  inner  layer  produces 
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For  the  third-order  outer-layer  slailarity  foraulation,  with  x,  r>  as  the 
Independent  variables,  the  dependent  variables  are 


*3(x,y)  -  *3(x,i?)  -  f3(h): 
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t3(x,y)  *  r3(x,jj)  **  [(l-2p*)~1/2  ]  $3(u)  . 


(6.19b) 
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With  p3(x,y)  -  P3(x.rj)  -  [(l-2pQ)  (dpQ/dx)  ]*3(r))  -  0. 
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The  following  boundary-value  problew  results: 
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Here,  It  Is  suggested  that  $  ~  i  i]  +  ...  as  t)  ■*  0. 

1  13 


To  suaaarize  the  results  obtained  for  the  outer-layer  siailarlty 
formulations,  it  is  convenient  to  Introduce  the  quantity 


a,  «  8  [(1-2PQ)'1  <^)  ] 


*  -1  dpo  2/3 

o'  <dT> 

=  6  [(Hq  ♦  2)  uJ0(x-x0)(H0+1)/H0+2)]~2/3, 


(6.23) 


such  that  ~  0(6)  «  1  for  x  fixed.  In  turn,  the  outer-layer  streaawise- 
velocity  and  Reynolds-stress  functions  take  the  following  forms: 
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Recall  that  un  -  (l-2p*)1/2  *  u*n(x  -  xrt)"1/(H0+2)  .  The  inner-edge  (tj  -  0) 


0  '  "r0'  “00'"  "O' 

behaviors  for  these  velocity  and  stress  functions  are 
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The  inner-edge  (jj  -*  0)  behaviors  for  the  velocity  and  stress  functions  are 
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With  the  Introduction  of  a4,  the  outer-layer  coordinate,  rj ,  and  the  inner- 
layer  coordinate,  C,  are  related  by 


7)  »  and/or  ?  -  a^  jj  . 


For  consideration  of  the  corresponding  inner-layer  velocity  and  stress 
functions,  it  is  convenient  to  introduce 
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Froa  (6.25)-(6.27) ,  it  is  determined  that  the  outer-edge  (?  ■*  «•)  behaviors  of 
the  inner-layer  velocity  and  stress  functions  are 
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6.2  Higher-Order  Approximation*  for  the  Inner  Layer 


The  equation*  of  motion  for  the  first-,  second-,  and  third-order 
approximations  for  the  inner  layer  are  presented  in  (4.32)  to  (4.34). 


For  the  first-order  approximation,  the  first  Integrals  of  the  momentum 
equations,  when  the  surface  boundary  conditions  are  taken  into  account,  are. 


Qj  -  fnc(x) , 
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From  matching  with  the  outer  layer  [see  (6.28)],  it  follows  that 
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Thus,  the  first-order  inner-layer  solutions  are 
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For  the  second-order  approximation,  the  first  Integrals  of  the  momentum 
equations  are 
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Matching  with  the  outer  layer  produces 
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*,<x.r)  -  0,(x.C)  -  [ (l-2pft)  (~)  ]  gJZ): 


•  -1  dpO  d*2 
*2  '  I<1_2pO)  dx  1  d{  * 


(6.34a) 


*  -i  dpo  4/3 


tj,(x, r)  -  t,(x.?)  -  ((l-2pn)  ]•*(?)  . 


(6.34b) 


*  dp  2/3 

q2(x.r)  -  q2(x)  -  -  »2  (— ) 


(6.34c) 


+  •  *  -4/3  0 

With  «2  »  [ (l-2pQ) (dpQ/dx)  JOa2/3r)  =*  const.,  the  second-order  inner-layer 


boundary-value  problea  is 


“2  ‘I  (H0  +  4)tV 


dg0  2 

d?l]  *  “2 


d«2 

-  0,  ®2  -  0  as  ?  -  0. 


2(V4)  ..3/2 

9«3  ^  ’ 


2<V4)  2 

- ; — C  as  (  -  »  . 

3kZ 


(6.35) 

(6.36a) 

(6.36b) 
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The  outer-edge  (C  -»  «•)  behaviors  for  the  velocity  and  stress  functions  are 


dg2  -  u+  ,3/2  .  +  _  .+  „l/2 

dT  b23C  +  b22<  *  b21?  + 


«  •  » 


2(H  +4) 

With  bM  '  4 - 3“’  b22  "  *12 


9* 


4 r  (hqm) 
9k2 


+2 


b21  "  *21  "  3*  (W2  +  2^2  1 2  )(H0+4) . 


+  2  +■  3/2 

*2  ’  *24<  +  *23C  +  *22C  +  ’ 


*  ♦  2<V4>  ♦  > 

Wlth  *24  "  *04  "  '  “  2  ’  *23  “  ♦is 


3* 


8 r  (Hq+4) 
9* 


+2 


*22  “  *22  “  3(W2  +  2  ^  )(H0  +  4) . 

K 


With  Bj,  n1  -  0,  the  first  Integrals  of  the  aoaentua  equations  for 


third-order  approximation  are 


q3  »  q3  »  fnc(x) , 


^3  ,fv° 

l3  dr  dx  r  *  3r 


fnc(x) 


Proa  matching  with  the  outer  layer,  it  Is  seen  that 


m3  -  0,  t3  -  0,  q3  -  q3  -*  0  as  r  -»  «. 


Thus,  the  third-order  inner-layer  solutions  are 


»3(x.r)  =  0,  t3(x.r)  -  0.  q3(x,r)  -  q3(x)  -  0. 


(6.37) 

the 

(6.38a) 

(6.38b) 

(6.39) 

6.40a-c) 


In  suanary,  then,  the  results  for  the  streaawlse-veloclty  and  Reynolds- 
stress  functions,  froa  the  Inner-layer  similarity  foraulations ,  can  be  expressed 


(6.40a) 


i  /v  _ w 

"  ~  -  1  x  A  ■  4.  • 

*  d?  *  d?  • ‘ •  * 


po 


•2 


=  +  . . .  . 


PO 


The  outer-edge  ({-►«•)  behaviors  for  these  velocity  and  stress  functions 


+  „l/2  .  +  >  -1/2  +  -l 

[b01?  +  b00  +  a01  <  +  a02C  + 


•] 


po 


.  ri  +  „3/2  . +  . +  1/2  , 

♦  Mb23C  +  +  b>iC  *  -1  + 


22  21 
+  „-l/2 


•2 


~  t  +  r  ^  *  -ift*  »-l/2  , 

l"02C  +  *00  +  +  -1 


or 


po 


,  +  2  +  3/2  + 

*  A>[*2«'  *  *23'  *  *22'  *  *  • 


In  (6.25)  and  (6.41). 
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02  02 
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.*  * 
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v04  24 


2(H  +4) 

0  .+  ♦ 

2  ’  *13  '  *23 


8r  (Hq+4) 


3k 


9k 


f(*L  ♦  )  (H„  +  4) . 


+2 


-  a 


32 


02  2*2 


*21  ”  b21  “  3*  <W2  +  2  V)(H0  +  4> . 

2c 


(6.40b) 


are 


(6.4la) 


(6.41b) 


(6.42a) 


(6.42b) 


• «  |  •  <  «  A  •  4 • » J  4  |  »  i  a. I ' ■  |  a  4 ’ a  % ' j 


la*  la»_L'  *.»  >.»  l.«  L*  jJV  W»  *->».*  |a*WV 


7.  SOLUTIONS  FOR  THE  DISTINGUISHED  INTERMEDIATE  LAYER 

In  Section  6,  for  the  streamwlse-veloclty  and  Reynolds-strees  functions, 
the  inner-edge  behaviors  for  the  outer-layer  solutions  and  the  outer-edge 
behaviors  for  the  inner-layer  solutions  are  given.  An  exaalnatlon  of  (6.25)  and 
(6.41)  suggests  that  a  distinguished  intermediate  layer  exists.  For  the 
distinguished  similarity  intermediate  layer,  the  independent-variable 
transformation  is  (x.y)  -  (x,x).  where 

,  _ .  1/2  -1/2  tl/2_ 

X  -  (»)C)  -  A,  I)  -  a,  C 


[(l-2p*)(^®)‘4/3J'I/2  <4"2r> 


(7.1) 


The  appropriate  speed  function  for  this  Intermediate  layer  is 


-lo  ■  <VWI/2  ■  »i/4  <  -  %1/4  V 


-  'I/4i'>-%»t^2»2/3)l/* 


21/4  V,  ■<«»  .  2,  -  V‘V2>"V2’ 


(7.2)  £ 


Introduction  of  z  •  A,  r\  into 


*  1/4  *  *2 

(u/u  )  -  A,  4  (U/U.n).  <r/u  Z) 


*1/2.  .  *2. 
A*  <r/uin)- 


and/or  x  «  a*^2{;  into 


•  -1/4  *  *2  -1/2  *2 

(u/up0)  *  A*  (u/ui0)>  (r/upo)  ■  A*  (r/ui0)* 


a*  given  in  (6.25)  and/or  (6.41),  produces 


u 


ioJ 


*10 


1/2. _+  3/2  +  -1/2.  .3/4..+  +  -1. 

♦  A*  IV,*  +  «21*  1  +  A,  I/S12*  ♦  «„X  ]  + 


12*  32 


(7.3a 


t  ^  ,  .  1/2.  ♦  2  ^+  .  3/4.  .+  3/2  +  -1/2  . 

"i?  ^  L+02X]  *  A*  [*04*  *  ^20  *  A*  [*13*  +  ^31*  J  + 

U10 


(7.3b 


The  results  of  Kader  A  Yagloa  (1978)  can  be  compared  with  those  of  the 
present  paper.  The  notation  of  Kader  A  Yagloa  In  terns  of  the  present  notation 
Is  jlven  by 


~2  m  * 

,  ,  ,1  dP,  ^0 

"  KY  “  (p  dx  KY  (  ~  1  dx  ’ 

c 


*  ~  •  A.-1 


(i,KY  -  c,ho  ■ ca  ^(1-2po)(i^,  1  * 


.  ,  . rw. 1/2  ~  .  +l/2.dp0.1/3 

“*  KY  *  p~  KY  •*  U-*»  <£T> 


(5?>KY  "  ^+Kl-2pn)  (~)  J  -  -  0; 


•  _j  %  2/3 


'dx 


1/2 


dp 


<y>KY  *  c  y:  (oy)^  -  u^  (—  yP  "  -  u_  u4„x 


0  1/2  ~  •  1/2 
U.  .1 
»  10 


«"»KV  *  “ 


(7.1 


In  turn,  the  Kader-Yagloa  representation  for  the  streaawlse  velocity  nay  be 
expressed  as 


[4.5  {1  ♦  10(«+A,)}1/2Z1/2] 


[U  »  )C{«y)1 


♦  a1/4[(3«+)1/2{2.44  log  r  -is  r  1/2  -  er  1u , 


where  r 


«n-vVi3/i 


5(1  *  10|»  *,)) 


The  results  of  Afzal  (1983)  also  can  be  compared  to  the  present  ones, 
notation  of  Afzal  in  terms  of  the  present  notation  is 


,  /V73  -  V'W73  -  ,2/3  • 

V*  (dT»  •  “J> 


•  d*>0  ■* 

\  / _ y  k  i 


(d)A  -  C5h0  -  cd  [(l-2p0)  (^)  ]  : 


<Va  -  ^A  -  *~4/3  f€i-.p:>-1«2*)a/ai-»  -  «-"V  . .. 


O'  dx 


rw  ♦  •  -1  dP0  2/3 

(A)A  “  (p~d^A  *  t(1"2P0)  (5T>  *  -  *  A.  *  0. 


(AR  )  -  d_1/V  - 

P  A 


(V)A  -  cy:  *  d'7/6  (£,l/V'«  - 


(u)  -  u  u:  (- — )  -  d"2/3A~1/4  - — 

*  'a  •  1 u  'a  *  * 

P  ui0 


In  Afzal 's  analysis,  it  is  taken  that  (aj  -  o.  (R  )  - 

2  A  p  A 

Note  that  (u,/ad)  -  (a)  -  •  a.  *  o. 


•.  and  (ARp) 


Afzal's  inner-layer  representation  for  the  streaawise  velocity  and  its 
expression  in  the  present  lnteraedlate-layer  notation  are 


Hp*  3.5  (1  +  5.4A) 

P 


U  y  1/2 

♦  2.5  (ARp)]A 


[3.5  {1  ♦  5.4  («+A,)}x1/2]  ♦  aJ/4[2.5(51/V)] 


(7.7) 


The  representation  of  Stratford  (1959b)  for  the  streaawise  velocity  in  the 
interaedlate  layer  is 

(“  ■  f  <;  ;j||>1/2y1/2  *  c  (*  £)1/31S.  -lth  C  -  const . 


-f--  tf  xI/2l  <•  *i/4tc  o2/3l 

ui0 


(7.8) 


Thus  the  lnteraedlate-layer  streaawise-velocity  representations  of  Kader  & 
Yagloa,  Afzal .  and  Stratford,  i.e.,  (7.5),  (7.7),  and  (7.8),  respectively,  are 
coapared  to  the  present  one.  (7.3a).  The  coaparlson  Indicates  that  the  three- 
layer  theory  provides  an  iaproved  characterization  or  the  overlap  doaaln  over 
that  provided  by  the  (classical)  two-layer  theory. 
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EXTERIOR  REGION:  Y-v  0(1); 

U  ~  0  (1),  T  ~  0 


DISTINGUISHED 

INTERMEDIATE 

LAYER:  y(  - Y/d372  -  y/d1/2  =* d1/2 r ~ O (1); 
U  «d1/4Uj,  T-d1/2r. 


INNER  LAYER:  r «  Y  /  d^  -  y  /  d  ~  0  (1); 
U-d1/2m,  T-dt 


Figure  1.  Schematic  diagram  of  the  asymptotic  structure  of 
the  boundary  layer. 


